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We present an ofF-resonant excitation scheme that reaUzes pronounced stationary inversion in a 
two level system. The created inversion exploits a cavity- assisted two photon resonance to enhance 
the multi-photon regime of nonlinear cavity QED and survives even in a semiconductor environment, 
where the cavity decay rate is comparable to the cavity-dot coupling rate. Exciton populations of 
greater than 0.75 are obtained in the presence of realistic decay and pure dephasing. Quantum 
trajectory simulations and quantum master equation calculations help elucidate the underlying 
physics and delineate the limitations of a simplified rate equation model. Experimental signatures 
of inversion and multi-photon cavity QED are predicted in the fluorescence intensity and second- 
order correlation function measured as a function of drive power. 
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The term "population inversion" refers to a system 
in which an excited state population exceeds that of a 
lower level. Its use is of broad interest to laser science. 
Semiclassical arguments readily show that a steady state 
drive, e.g., a continuous wave (cw) laser, cannot cre- 
ate population inversion in a two level system (TLS) 
because of stimulated emission; thus, standard atomic 
and molecular lasers typically exploit a cascaded through 
three, four or more levels to achieve inversion In 
1988, an important limitation of the semiclassical argu- 
ment was pointed out by Savage [2], who demonstrated 
that the multi-quanta regime of an atom-based cavity 
QED system can achieve stationary population inver- 
sion in a TLS. Specifically, by including the effects of 
quantum noise^ and with the pump resonant with the 
degenerate atom-cavity resonance. Savage [2] — and later 
Lindberg and Savage [3] — predicted achievable inversions 
of around 0.01-0.03. The unanticipated result was ex- 
plained in terms of a 4-state two-quanta model of an on- 
resonance cavity-atom system in the bad-to-intermediate 
cavity regime. Lindberg and Savage (sj suggest that the 
inversion might be measured in high pump resonance flu- 
orescence, as a small excess (a few %) over the saturation 
fluorescence intensity. To the best of our knowledge this 
prediction has not been verified by experiment on either 
atomic or semiconductor systems. 

There have been numerous recent experiments in semi- 
conductor cavity systems using coherent excitation, most 
exploiting a deliberate detuning between the cavity and 
quantum dote (QD). For example, Majumdar et al 01 
studied the coupling between a photonic crystal cav- 
ity and an off-resonant QD under coherent excitation of 
either the cavity or the QD; they investigated exciton 
power broadening under coherent excitation in the pres- 
ence of a detuned cavity. Ulhaq et al [HI] investigated 
linewidth broadening of a resonantly excited QD in a 
micropillar cavity, again with QD-cavity detuning. Both 
groups demonstrate significant power broadening of the 



target QD exciton, and, importantly, the experimental 
capacity to explore off-resonant QD-cavity systems using 
nonlinear coherent excitation in a controlled way. There 
have also been remarkable developments in circuit QED 
0-[9|, where anharmonic cavity QED effects are realized 
at dipole coupling strengths two to three orders of mag- 
nitude larger than in atomic or QD systems. 

In this Letter, we introduce a scheme to create pro- 
nounced population inversion in a dissipative TLS cou- 
pled to a off-resonant cavit y. Supp orted by quantum tra- 
jectory (QT) simulations [10|, [uj and quantum master 
calculations, we show that off-resonant excitation of a 
detuned TLS-cavity system, with intermediate to strong 
couphng, results in population inversions that can eas- 
ily exceed 0.25. These stationary inversions are enabled 
through two-photon excitation of the dressed states, and, 
strategically, only occur for a detuned TLS and cavity. In 
fact, significant prototype inversions have been seen in a 
circuit QED device excited via so-called sideband transi- 
tions [l2|, an excitation scheme similar to ours; however, 
the unanticipated inversion is discussed only briefly in 
terms of a simplified rate equation model. We use QT 
simulations to, first, connect to, then extrapolate beyond 
rate equations. We show that they have a limited range 
of validity, and focus on the possibility of achieving sig- 
nificant inversion far outside the rate equation regime. 
We show that such inversion exists for practical semicon- 
ductor systems where it is potentially important for mi- 
crolasers and few photon quantum light sources. Other 
non-standard routes to gain include lasing without in- 
version [l3|, but here we are dealing with conventional 
inversion-driven gain, but with an unfamiliar pumping 
scheme. We also suggest experiments for measuring the 
inversion, via the fluorescence intensity measured as a 
function of drive power 0, [H, or through the second- 
order correlation function; both exhibit clear signatures 
of the proposed off-resonant stationary inversion regime. 

The cavity QED system is described through a quan- 
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FIG. 1: (Color online) (a) Quasienergy levels truncated to 5 
states (S = 6cx = coc — oox)- |1, +) refers to one created cavity 
photon with one exciton (excited TLS) created; for a far- 
detuned cavity and exciton, this doubly-excited state relates 
to the Jaynes-Cummings dressed states through |2L) ^ |1, +)• 
Bold arrows illustrate two-photon resonance, with the pump 
detuned from the cavity by about —6/2. Small arrows show 
the direction of the Stark shifts, (b) Simplified rate equation 
model, where Vp = Qtp/2K is an effective two-photon pump 
rate, with Qtp the two-photon Rabi frequency (see text). 



turn master equation for the reduced density matrix of 
the TLS (exciton in a semiconductor dot) and cavity. 
Working in the interaction picture, i.e., in a frame rotat- 
ing at the pump frequency, ujl^ the master equation is 

M 



dp 
~dt 



n 



[Hs,p] + ( KC[a] + lc[cj_] + ^C[cj,] ] p, (1) 



7' 



with system Hamiltonian 

+ fl7?^(cr++cr-), (2) 

where r]x is the coherent pump rate of the TLS, g is 
the TLS-cavity coupling strength, a is the cavity mode 
annihilation operator, <j+, <j~ and are Pauli operators, 
and SaL = ^a—^L are pump laser detunings; dissipation 
is described by the Lindblad C[£]p = '^^p^^ ~ ^^^P~ P^^^i 
with separate terms to account for decay of the cavity 
mode (rate a^), and radiative decay (rate 7) and pure 
dephasing (rate 7') of the TLS. 

Figure [11(a) depicts the dressed quasienergy levels of 
Hamiltonian ([2]), truncated to 5 states. The lowest 4 
levels are the main ones responsible for cavity-assisted 
inversion, where decay of |1, +) populates |0, +) through 
the emission of a cavity photon. Importantly, two photon 
resonance to |1,+), with negligible excitation of |1,— ), 
is achieved by setting 26 lx ^ ^ = ^cx = — ^x ^ 
5f,/^,7,7^ For increasing pump strengths. Stark shifts 
lower the resonance frequency. For comparison, a simpli- 
fied "rate equation" model, with adiabatic elimination of 
intermediate states, is shown in Fig. [T](b). 
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FIG. 2: (Color online) Sample QTs near, (a) [g = IOOOac, 
r]x = IOOk, 7 = 0.05k], and outside, (b) [g = 100/€, r]x = 150k,, 
7 = 0.05k], the rate equation regime. The blue solid curve 
plots the expectation of the TLS excited population (exciton 
population) conditioned upon the record of quantum jumps, 
while the red solid curve is the conditioned expectation of 
the cavity photon number. Quantum jumps are indicated by 
the crosses and circles (and vertical dashed lines) for photon 
emission via the cavity and TLS, respectively. The steady- 
state expectation of the exciton population is shown with the 
horizontal solid line, while the rate equation estimate is shown 
with the dashed horizontal line. The insets enlarge the regions 
near tK ^ 46 — 48 and tK ^ 29 — 31 in (a) and (b), respectively. 



We first derive the TLS population within the simpli- 
fied model. With coherent drive, and neglecting pure 
dephasing, we obtain the steady-state population of the 
TLS excited state (population of excitons), 
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with two-photon drive Qfp = {^Vx I / {d / • Note 
that this expression differs from the one used by Leek et 
al. [12]; their formula, fix = Vttp/ {^tp^l /'^f^^l^tp/'^f^)^ 
is obtained by treating Vttp as the rate of an incoherent 
drive. 

Using QT simulations we identify a regime where the 
rate equation picture is approximately valid while also 
allowing for significant inversions. Broadly, we require: 
6/n > 10^, g/n > 10^ and rjx/hi > 10^ Figure [2]; a) 
shows a sample QT on the boundary of this region, with 
j/k, = 0.05. Taking level shifts into account, the two- 
photon resonance is near 6xl = —0.49^. We highlight 
the following: (i) the TLS jumps to an excitation close 
to unity with the cavity excitation close to zero; (n) each 
such jump is initiated by the emission of a photon by the 
cavity corresponding to the transition |1,+) |0,+) 
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in Fig. [IJa); {in) all excitation periods terminate with 
the emission of a photon by the TLS corresponding to 
the transition |0,+) \^^—) in Fig. [TJa); (iv) the in- 
tervals between excitation periods (e.g., t/^ ^ 15 — 22) 
see the cavity and TLS populations grow continuously 
and in step, corresponding to coherent excitation of the 
|0, — ) transition. Solid and dashed horizontal 

lines show that there is nevertheless a discrepancy be- 
tween the exact (solving the ME in a suitably large pho- 
ton number basis 0) population {fix ^ 0.77) and that 
predicted by Eq. ([3]) {fix ~ 0.84); most likely this is due 
to our neglect of level shifts when deriving Qtp- 

We next choose parameters that start to deviate away 
from the rate equation picture. Figure [2fb) is plotted 
for S/k = 10^, g/hz = 10^ and rjx/i^ = 1.5 x 10^ stih 
with = 0.05. This places us close to the circuit QED 
regime of Leek et al. [l2|. Although a similar jump pat- 
tern is seen, the following differences are noticed: {i) the 
jump "up" reaches an excitation probability significantly 
less than unity (a small piece of |1, t/) ~ \^-,~) is mixed 
into the state realized after the jump); {ii) the drive 
strength has been increased, which induces Rabi oscilla- 
tions on the |0, — ) |1,+) transition (e.g., right before 
tn = 31); {in) the TLS and cavity do not continue indefi- 
nitely to Rabi oscillate in step, as for ttz ^ 26 — 26, where 
the TLS shows a damped evolution toward a number near 
0.9, while the cavity photon number decays to zero. This 
behavior illustrates a marked difference between QTs and 
rate-equation-inspired quantum jumps. The final state of 
the damped Rabi oscillation is \l^L) ^ |0,+); but this 
state is reached via a coherent evolution in response to 
the back action of a null measurement — the system lo- 
calizes to this state because no photon happened to be 
scattered through the cavity; the two-photon resonance 
still assists passage to the inverted state \l^L) ^ |0,+), 
but through the absence rather than the presence of a 
quantum jump. 

QT simulations have clarified two keys points: (a) the 
rate equation model is extremely restricted even for cir- 
cuit QED parameters, and (b) significant population in- 
version can be achieved when dissipation rates increase 
relative to the dipole coupling. This suggests that signif- 
icant inversions might be achieved using semiconductor 
QD parameters. 

Turning now to semiconductor QED systems, we adopt 
the following scaled parameters in the strong coupling 
regime: = ^/2.5, 7 = ^/200, and 7' = ^/25. For cur- 
rent strongly coupled QD systems, g is around 0.1 meV 
{^ 24 GHz), with similar values measured in a wide 
range of devices, including exaniples comprised of pho- 
tonic crystals [16| , micropillars [17| , and microdisk cav- 
ities [18|. To connect with measurements of intensity 
linewidth broadening (e.g., Refs. 0, Q), we calculate the 
intensity of the exciton mode, Ix cx n^^, or the cavity 
mode, Ic (X Uc. We solve the full ME, Eq. ([1]), in a basis 
that can be truncated at an arbitrary photon/exciton 
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FIG. 3: (Color online) Semiconductor quantum-dot QED sys- 
tem, (a-b) Pump-dependent steady-state exciton number ver- 
sus detuning calculated with 1-photon and 3-photon trunca- 
tions; for Tjx = 0.02g — 2g (blue to purple) and S = lOg. 
Note the clear peak in the 3-photon case near lul — (^x ^ S/2. 
(c) Sample QT simulation on the peak labeled pi; symbols 
are defined as in Fig. 2, with the additional pure dephasing 
phase-glitch indicated by an inverted triangle. 



state, allowing us to compute both weak-excitation- 
approximation results (truncation at one quantum) and 
the regime of multi-quanta cavity QED. While weak ex- 
citation is commonly assumed when analyzing dissipa- 
tive semiconductor cavity systems, we will show that the 
approximation drastically breaks down in the proposed 
excitation regime. 

We display in Figs. [3fa-b) the exciton population as a 
function of laser-exciton detuning, with S = 10 g^ for two 
different truncations. All calculations predict substantial 
power broadening and show the expected peaks at the 
bare cavity and exciton resonances, in agreement with 
recent weak-coupling experiments 0, With multi- 
photon effects included, however, and for sufficiently 
large drive, the two-photon resonance emerges approxi- 
mately half way between Ux and cjc- The resonance yields 
pronounced exciton populations, exceeding Ux = 0.75 
for r]x ^ 2g (inversions of w = Ux — 0.5 > 0.25). The 
position of this intermediate "resonance" changes with 
the pump strength, which gives a direct signature of the 
Stark-shifted two photon resonance. Similar peaks ap- 
pear in the cavity intensity (not shown). As a strong 
contrast, on-resonance excitation {S = 0) yields only very 
small inversions (a few %) UM- 

What is surprising about these predictions is that they 
show a pronounced influence from higher lying ladder 
states, which are notoriously difficult to see in a semi- 
conductor system because of the relatively large cavity 
broadening. Recent experiments exploring QD anhar- 
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FIG. 4: (Color online) (a) Exciton (blue) and cavity (red) 
population, and (b) second-order correlation function versus 
drive strength, r/x, for the semiconductor case of Fig. [S] 



monicities have been reported by Kasprzak et al [Toj , 
using advanced coherent nonlinear spectroscopy. The sig- 
natures observed are very weak. On the other hand, the 
scenario we investigate shows a profound impact from the 
second rung of the Jaynes-Cummings ladder, specifically 
the state |2— ) ^ \'^^-^)^ where a quantized level of the 
radiation field serves as a rung in a multi-level pumping 
scheme. Indeed, we stress that the two photon resonance 
and resulting inversion would not exist except for the in- 
fluence of the higher lying Jaymes-Cummings levels. This 
is an important flnding for assessing the nonlinear cavity 
regime. Even in atomic cavity QED systems, signatures 
of the two-photon resonance with resonant atom-cavity 
coupling are weak [20j. 

Having demonstrated how to achieve large population 
inversion through coherent off-resonance excitation in a 
dissipative semiconductor system, we now consider some 
common experimental techniques that might detect it. 
Savage [2] suggests that a maximum inversion created 
by driving at saturation would yield a direct signature 
in fluorescence. A method more speciflcally sensitive 
to two-photon resonance involves carrying out a photon 
cross-correlation, i.e., measuring the correlation function 
^^^c = {(J~^ cr~ a) = which allows one to access 
information about the quantum statistics and penetrate 
the pumping cycle that creates inversion. Figures Hl^a-b) 
display the dependence of the populations and second- 
order correlation function on the coherent drive strength 
for a semiconductor system. Two regions of maximum in- 
version are seen in the populations, with maxima in the 
cross-correlation function occurring in one to one corre- 
spondence. 

Finally, it is worth emphasizing that the predicted in- 
versions survive over a very wide range of experimental 
parameters — e.g., for the semiconductor system, when 7, 
7', and S all increase by a factor of flve, we still obtain 
peak populations ranging from 0.65 — 0.85. In fact, larger 
detunings can result in larger inversions since this avoids 
exciting the intermediate states, though, of course, more 
pump power is required. A more detailed coverage of 
the parameter dependence, including the influence from 
electron-phonon coupling with nonlinear coherent driv- 



ing 111, [iH, will be described elsewhere. 

In conclusion, we have proposed a cavity QED scheme 
with off-resonance excitation that achieves stationary 
population inversions exceeding 0.25 in semiconductor 
QDs. We also explored the connection to previous mea- 
surements in circuit QED and uncovered the limitations 
of a simplifled rate equation model: quantum trajectories 
clarify the dynamic leading to quantum-noise-induced in- 
version for dipole couplings ranging over three orders 
of magnitude in units of cavity loss. More generally, 
we demonstrated the existence of a cavity-assisted two- 
photon resonance that signiflcantly enhances the regime 
of multi-photon cavity QED, in particular for QDs. This 
can beneflt the continued worldwide search for genuine 
quantum signatures of cavity QED in a number of inter- 
esting material systems, along with the development of 
semiconductor-based single QD lasers (isj . 
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